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1. Introduction 



Abstract. EPOXI (EPOCh + DIXI) is a NASA Discovery Program Mission of Opportunity 
using the Deep Impact flyby spacecraft. The EPOCh (Extrasolar Planet Observation and Char- 
acterization) Science Investigation will gather photometric time series of known transiting exo- 
planet systems from January through August 2008. Here we describe the steps in the photometric 
extraction of the time series and present preliminary results of the first four EPOCh targets. 
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The high-precision time series produced from the EPOCh observations, with their sub- 
stantial time and phase coverage, allow us to pursue three important science goals with 
£> respect to transiting exoplanets. First, we will search for evidence of reflected light from 

k> ' the giant transiting planets at secondary eclipse (Rowe et al. and Matthews et al., these 
proceedings). Second, EPOCh will be sensitive to terrestrial planets in these systems, 
using two detection methods. The first will be a search for direct transit events caused 
by approximately Earth-sized planets, particularly "hot Earths" that are predicted to 
be captured in low order mean motion resonances with the giant planets (Croll et al. 
2007a, b). EPOCh will also perform transit timing measurements for the giant planets, 
and search for timing variations due to terrestrial planet perturbers (Agol et al. 2005; 
Holman & Murray 2005). The third science goal of the EPOCh mission is a precise 
determination of the system parameters, notably the planet radius, and a search for 
planetary satellites and circumplanetary rings through their perturbations to the transit 
light curve (Brown et al., 2001; Barnes & Fortney, 2004). These proceedings will sum- 
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marize the steps we take to produce the time series, as well as describe the preliminary 
search for additional transiting bodies. 

Each time series consists of continuous 50-sccond integrations through a clear visible 
filter with the 30 cm aperture High Resolution Imager (HRI) instrument. These obser- 
vations have a typical duration of 3 weeks per target. The spacecraft will then continue 
en route to its 2010 encounter with comet Hartley-2, the target of the DIXI Science 
Investigation. For a history of the EPOCh observations and more information about the 
EPOXI mission, please see Christiansen et al. (these proceedings). 



2. Data Reduction 

We use the existing Deep Impact reduction pipeline to perform bias and dark sub- 
traction, as well as flat fielding (Klaasen et al. 2005). We then conduct the photometric 
extraction as follows. We determine the position of the star on the CCD using PSF fit- 
ting, by maximizing the goodness-of-fit (with the x 2 statistic as an estimator) between 
an image and a model PSF with variable position, additive sky background, and mul- 
tiplicative brightness scale factor. The PSF itself, with a roughly 10 pixel FWHM, is 
produced from the drizzle of more than 1200 60x60 pixel cutouts, filtered to eliminate 
cosmic ray hits before drizzle. The final PSF product is sampled to a tenth of a pixel. 
A bilinear interpolation of this PSF increases the sampling to a hundredth of a pixel, 
which is the accuracy to which we estimate the position from the x 2 grid. At this point, 
we perform cosmic ray filtering by removing images from the sample with a larger than 
lOcr outlier in the residuals between image and best-fit PSF model. Because of the high 
cadence of the EPOCh observations, we simply reject any images containing a cosmic 
ray overlying the stellar PSF from the time scries. 

Flat fields for the High Resolution Imager were last taken while the spacecraft was 
on the ground. We construct secondary flat fields using a "stimulator" LED on board. 
Because this LED is mounted to one side of the CCD, we first remove an overall trend 
of increasing brightness with position by dividing out a best-fit plane to the flat field. 
We then produce a final "stim" by coadding the 100 stims taken in each block of stim 
observations, using pixel-by-pixel outlier rejection to account for cosmic ray hits. We 
observed changes to the stim in the form of dark spots continuously during observations; 
we attribute this to possible radiation damage. However, the standard deviation of the 
time series was unchanged when we removed observations for which one of these dark 
spots was present in the aperture. 

After we use PSF fitting to determine the position of the star on the CCD, we use 
aperture photometry to perform the photometric extraction. The structure of the CCD 
includes two rows in the middle of the array with physically smaller pixels, such that the 
PSF has a different shape when the star straddles these rows. We account for this problem 
by scaling these rows by the multiplicative factor that best recovers a flat dependence of 
out-of-transit flux on row position. We also perform a similar multiplication by a scale 
factor for the two middle columns. In the first version of the photometry pipeline, we 
also divided out a spline fit to the flux, first as a function of y position, and then as a 
function of x, to account for additional flux dependence on position. We decided on an 
optimal aperture radius based on analysis of the standard deviation of the out-of-transit 
time series. We found that this standard deviation was minimized for an aperture radius 
of 10 pixels, corresponding to twice the HWHM of the PSF. 

After performing the stim division and the scaling of the two middle columns and rows, 
the time series still suffered from significant red noise. At this point, we implemented a 
2D spline fit to the data with position, which largely removed the red noise when divided 
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Figure 1. Top panel: HAT-P-4 time series before (lower curve) and after (upper curve) 2D 
spline correction. The uncorrected time series has had the two middle rows and columns in 
each image scaled by a multiplicative factor to reduce the flux dependence on position, but the 
images were not divided by a flat field constructed from a stim. Instead, we have used the 2D 
spline to correct for interpixel variations. Bottom panel: The data (diamond symbols) bin down 
consistently with the expectation for Gaussian noise (shown with a line, renormalized to match 
the value at N=l). 

from the time series. Figure 1 shows the time series before and after this 2D spline 
correction. In the bottom panel, we show that the time series after the 2D spline bins 
down roughly as predicted for Gaussian noise. In the corrected time series, we attribute 
the scatter to photon noise and low-level cosmic rays. 



3. Preliminary Results and Discussion 

We have observed seven transits of HAT-P-4, five transits of TrES-3, three partial 
transits of XO-2, and seven transits of GJ 436. Figure 2 shows the EPOCh observations, 
phase-folded and binned, for these four systems. Our strategy is to conduct the science 
analysis of the data while working in parallel to refine the photometric precision through 
improvements to our extraction pipeline. We are currently conducting a Markov Chain 
Monte Carlo modeling and reflected light analysis to refine the system parameters and 
constrain the planetary albedos. These analyses are explained in more detail in Chris- 
tiansen et al. (these proceedings). We are also conducting preliminary searches for shallow 
transits made by additional bodies. For GJ 436, this search is especially interesting, since 
the eccentricity of the known transiting planet points to the possible presence of a second 
planetary companion. We are just beginning to develop the software to search for these 
additional transits, which will take steps briefly outlined here. 

The set of lightcurves for additional transiting bodies is a three parameter family if we 
make the simplifying assumption of no limb darkening of the host star and an inclination 
angle i of 90°: radius of the planet R p , period of the planet P, and phase (p. Using the 
Mandel and Agol (2002) routines for generating lightcurves given these parameters, we 
compute a grid of models corresponding to additional possible planets. For each of these 
models, we compute the x 2 , using the out-of-transit standard deviation to estimate the 
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Figure 2. Phased and binned light curves for the four EPOCh targets observed to date. We have 
applied the 2D position decorrelation to HAT-P-4 and TrES-3. The solid curve corresponds to a 
best-fit transit curve solution to these data, using Mandel & Agol (2002) routines for generating 
light curves. The limb darkening coefficients were assumed for HAT-P-4 and TrES-3 from stellar 
models, but were allowed to vary for XO-2 and GJ 436. 

error in each point. At the positions of local minima, we will change to a routine that 
more finely samples the % 2 space so that we can accurately gauge the significance levels of 
the parameters. The improvement in the x 2 over the null hypothesis will also determine 
whether we can claim a positive detection. The photometric precision of the GJ 436 time 
series (with S/N=44 for each transit event) would enable a detection of a 1.3 R ffi planet 
with S/N of 4 and a 2.1 R e with S/N of 10, even if the planet produced only a single 
transit event. 
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